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We study the co-evolution of Yang-Mills fields and perfect fluids in Bianchi type I universes. We 
investigate numerically the evolution of the universe and the Yang-Mills fields during the radiation 
and dust eras of a universe that is almost isotropic. The Yang-Mills field undergoes small amplitude 
chaotic oscillations, which are also displayed by the expansion scale factors of the universe. The 
results of the numerical simulations are interpreted analytically and compared with past studies of 
the cosmological evolution of magnetic fields in radiation and dust universes. We find that, whereas 
magnetic universes are strongly constrained by the microwave background anisotropy, Yang-Mills 
' universes are principally constrained by primordial nucleosynthesis and the bound is comparatively 

£\j , weak, and J7ym < 0.105f2 ra d . 
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I. INTRODUCTION 

. There is considerable interest in the generation and evolution of magnetic fields in cosmological models 0. This 
interest has a double focus. There is the hope that an explanation might be found for the existence of significant 
magnetic fields in galaxies. This may require seed fields to originate very early in the history of the universe, or even 
■ form part of the initial conditions. We have also discovered that the cosmological evolution of anisotropic stresses, 
like electric and magnetic fields, has interesting mathematical features that create unexpected physical consequences. 
These features appear when a homogeneous and anisotropic stress with a trace-free energy-momentum tensor is 
present during the radiation era of the universe. The isotropic black-body radiation and the anisotropic stresses 
evolve to first order in the same way, but there is a greatly slowed (logarithmic) decay of the shear anisotropy caused 
' by the pressure anisotropy of the anisotropic fluid [J, HI- If the evolution were to be linearised around the isotropic 
model then a zero eigenvalue would exist. When the equation of state of the accompanying perfect fluid changes from 
radiation to dust the evolution of the shear is still dominated by the pressure anisotropy of the anisotropic stresses 
{— j-i but the zero eigenvalue disappears and the shear falls as a reduced power of the cosmic time. In the study of these 
effects made in ref Q we considered the effects of anisotropic stresses whose trace-free anisotropic pressure tensor ix a b 
was proportional to a density /j,a, so 
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TTab = CabfJ-A, (1-1) 

with C a b constant and C° = 0. Barrow and Maartens considered the generalisation of these results to general in- 
homogeneous cosmological models close to isotropy using the covariant formalism [5( and also some applications to 
Kaluza-Klein cosmologies 0. The key feature in the evolution of almost-isotropic cosmological models containing 
perfect fluids plus a magnetic field is that during the radiation era the ratio of the shear to Hubble expansion rate falls 
only logarithmically in time. The universe evolves towards an attractor state in which it is proportional to the ratio 
of the magnetic and perfect-fluid densities. This means that the ratio of these densities at the epoch of last-scattering 
determines the large-angle temperature anisotropy of the cosmic microwave background (CMB). This enables us to 
use observations of the CMB to place strong bounds on any homogeneous cosmological magnetic field Q , or other 
anisotropic stress defined by a constant C a b, 4]. An important extension of these studies is to consider anisotropic 
stresses with more general pressure tensors, for example those with time-dependent C a b- 

An important case with time-dependent C a b, which generalises the electromagnetic field, is that of a Yang-Mills 
stress. This has been studied for pure Yang-Mills fields and no accompanying fluid by several authors in simple 
Bianchi t ype I anisotropic universes @,|UEi> m au class A Bianchi type universes ^J, and in Kantowski-Sachs 
universes [T^l • These studies reveal that the Yang-Mills field creates a form of chaotic behaviour during the evolution 
of the Yang-Mills stress |lfjj | . This is not surprising because such behaviour is present in Minkowski space-time and 
is nothing to do with the cosmological evolution or the spacetime curvature. It should not be confused with the 
chaotic behaviour of general relativistic origin that is found in the Mixmaster' universes of Bianchi types VIII and 
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IX, even in vacuum 1 . Therefore the Yang-Mills chaos occurs in the simplest axisymmetric universes of Bianchi type I 
and is present even when the anisotropy level is arbitrarily small. This raises the interesting question of whether the 
chaotic evolution might leave an imprint in the temperature anisotropy of the CMB when it decouples at a redshift 
z rec ps 1000. 

However, all of these earlier studies consider the evolution of Einstein- Yang-Mills equations only for the case of a 
pure Yang- Mills stress in an anisotropic cosmological model. Our experience Q with the behaviour of magnetic fields 
in anisotropic universes teaches us that it is important to include the presence of a perfect fluid in order to find the 
realistic evolution of the Yang-Mills stress during the radiation and dust-dominated eras of the cosmological expansion. 
In this paper we consider this generalisation by analysing the evolution of spatially homogeneous, anisotropic universes 
of Bianchi type I containing both a perfect fluid and an anisotropic Yang-Mills stress. By combining analytic and 
numerical studies we determine the evolution of the shear and the Yang-Mills stress during the early universe in order 
to discover if it is possible for it to leave an observable signature in the CMB. 

In section 2 we set out the Einstein- Yang-Mills equations for the Bianchi class A universes. In section 3 we 
specialise to consider in detail the evolution for the Bianchi type I universe containing Yang-Mills fields and perfect 
fluids, focussing on the physically significant cases of pressureless dust and black-body radiation. In section 4 we 
describe the numerical and analytical results and in section 5 compare them in detail with the situation that results 
when a pure magnetic field replaces the Yang-Mills field. We consider the observational bounds that can be placed 
on the magnetic and Yang-Mills fields in the presence of perfect fluids using the microwave background temperature 
anisotropy and the primordial nucleosynthesis of helium-4 in section 6. We summarise our results in section 7. 



II. THE EINSTEIN- YANG-MILLS EVOLUTION EQUATIONS 

The Einstein equations to be solved for universes containing a perfect fluid and a Yang-Mills (YM) field are 

Gfj, v — (2.1) 

T^u = Tyyiixu + T m)iu , (2-2) 

T mMV = diag(- /i m , Pm,Pm,Pm), (2.3) 
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F (A) p (A)\_\ F (A) F (A)\a 



(2.4) 



where FjuP is the field strength of YM field, 11 A" describes the components of the internal SU(2) space, and /i m and 
p m have an equation of state 

Pm = (7 - !)A*m- 

We shall assume that Tym^v and T mliv are separately conserved and there is no energy exchange between them. 
Hence, from the vanishing covariant divergence of Tym/^ we find that the YM equations are 

F^. v -ixF" = 0, (2.5) 

where c?ym is the self-coupling constant of YM field, and F^ lL/ — F^^ta and A v — A^ta with ta being the SU(2) 
basis. Setting our units as 8irG/g Y M = 1 2 an d c = K = 1, the basic equations become free of the value of <7ym- 

We adopt the orthonormal- frame formalism, which has been developed in ref. |18| . In Bianchi- type spacetimes, 
there exists a 3-dimensional homogeneous spacelike hypersurface E t which is parameterized by a time coordinate t. 
The timelike basis is given by eo = dt- The triad basis {e a } on the hypersurface £4 is defined by the commutation 
function "f c ab as 

[e a , e b ] =i c ab e c . (2.6) 



1 The most general cosmological analysis of the Yang-Mills chaos is that carried out for all Ellis-MacCallum Class A Bianchi-type vacuum 
cosmo logi es by Jin and Maeda 111! and they are able to consider the simultaneous presence of the Yang-Mills chaos and Mixmaster 

chaos Ea, m, EE EE E3i 

2 This choice of units has the following consequences. If <?ym = b then 8irG = 1; if gYM = 10 — 10 , then 8nG = 10 -20 and in this case the 
unit of time is 10 10 t p i ~ 10~ 34 sec and the present Hubble time is H^ 1 = 1.2 X 10 10 yr = 4 X 10 17 sec = 10 51 . 
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It is convenient to decompose 7 c ab into the geometric variables denoted conventionally by the Hubble expansion H, 
the acceleration u a , the shear a ab , and the vorticity uj ab , an d the variables f2 a (the rotation of e a ), and the variables 
a a and n abl which distinguish the type of Bianchi model, so 

lib = -v b a -H6 b a -e a bc (u c + n c ), (2.7) 

7°oa = u a , (2.8) 

7° a6 = -2e o6 c a, c , (2.9) 

l c ab = e abd n dc + a a S b c -a b S a c . (2.10) 

If Bianchi spacetimes are expressed in this way, the vorticity uj ab always vanishes because they are hypersurface 
orthogonal. 

In order to analyze the EYM system, we shall study the simplest case. First, using the gauge freedom of the YM 
field, we set Aq (t) = 0, which simplifies the vector potential to A = A a A \t)TA'^ a , where u> a is the dual basis of 
e a []Jj- Next we shall restrict attention to cosmologies of Bianchi type class A, in which the vector a a vanishes. We 
can also diagonalize n ab , i.e. n ab — diag{ni, ri 2 , n^}, using the remaining freedoms of a time-dependent rotation of 

the triad basis. Then, we can show that if a ab , Aa and A a A ^ do not have off-diagonal components initially, the 
equations of motion guarantee that those variables will remain diagonal during the subsequent evolution for the class 
A Bianchi spacetimes (see the Appendix in |ll|). As a result, f2 a vanishes and the number of basic variables defining 
the initial value problem reduces from 21 (12 for spacetime [H, N a , a ab , f2 ] and 9 for YM field [A^]) to 10 (7 for 
spacetime [H, N a , a aa ] and 3 for YM field [Y ])■ 

In order to discuss the dynamics most efficiently, it is convenient to introduce the Hubble- normalized variables, 
which are defined as follows: 

£ Qb ^, = (2.11) 

We also re-express the shear variables of the Bianchi spacetime as 

£+ = - (£ 22 + £ 33 ) , £_ = — (£22 ~ ^33) , S 2 = Y? + + Y?_ . 
The diagonal components of the YM field potential are described by new variables a(t), b(t), and c(t) as 

a^4 1} , b^A?\ c^Af. 

and a new time variable, t, defined by dr = Hdt is introduced so that t denotes the e-folding number of the scale 
length. 

Using these variables we can write down the evolution and constraint equations explicitly. They consist of the 
generalized Friedmann equation, the dynamical Einstein equations and the YM evolution equations. The generalized 
Friedmann equation, which is the constraint equation, is 

s 2 + n YM + n m + k = i, (2.12) 

where !)ym is the density parameter of the YM field, i.e. the Hubble-normalized energy density of the YM field, 
which is defined by 



1 



G 



{a' + (-2S+ + l)a} 2 + \b' + (S+ + V3S_ + + |c' + (E+ - \/3S_ + l)cX 



be \ 2 / , ca\ 2 ( ab^ 



bc\ t ca 

NiCL+ h) +( iV2fe + 77 



A 3 c+- 



(2.13) 



f2 m is the density parameter of perfect fluid, defined by f2 m = /^ m /(3_ff 2 ), and K is the Hubble normalized curvature, 
which is defined by 

(3)R 1 

K = = — {N 2 + A 2 + A3 2 - 2(AxA 2 + A 2 A 3 + A 3 7Vi)} . 
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Note that the positive 3-curvature corresponds to K < 0. Hence, except for the Bianchi type IX models we have 
K > 0. The energy density of YM field and perfect fluid are always positive definite. Thus we find that £ 2 , SIym, 
and K are restricted to the domains < £ 2 < 1, < SIym < 1, and < K < 1, except for the Type IX universes. 
The dynamical Einstein equations are 



H' 

K 
£'_ 

K 



( g -2)E + -5++n +! 

(?-2)E_ -SI + n_, 

(g-4E+)iV 1 , 

(<Z + 2£+ + 2\/3£_)7V 2j 



where q is deceleration parameter: 



2£ 2 



in. . 

2 2# 2 



= 2£ 2 



YM + -(3 7 -2)fi B 



The quantities S± are given only by the N a : 

1 



S+ 
S- 



G 



{(N 2 - N 3 ) 2 - JVi(2JVi — N 2 — N 3 )} 
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2^ 



(Na-NiXNx-Nz-Ns). 



(2.14) 
(2.15) 
(2.16) 
(2.17) 
(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



Last, we define the Hubble-normalized anisotropic pressures of YM field, II + and II_, so that 
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-2 {a' + (-253+ + l)a} 2 + {&' + (£+ + \/3E_ + l)b} 2 + {V + (E + - V3£_ + l)c}' 



-2 ( 7Via + — 



iV 2 6H 



iV 3 cH 



(2.23) 



n_ = — 



V3 
G 



{&' + (£+ + V3E_ + l)b} 2 + [c' + (£+ - V3£_ + l)cV 



iVaC- 



aft 



(2.24) 



Here we can see the complexity introduced by the YM dynamics. We can see how the anisotropic pressures are 
driven by the directional scale factors and their time derivatives, and not solely by the fractions of the density as was 
the case with the simple model of anisotropic stresses defined by eq. iJOJ. 

The isotropic matter conservation equation is 



^ = {2g-(3 7 -2)}0 B 



(2.25) 



and the YM evolution equations are 
a" = (q-2)a' 



4E+ -Nx{N 2 + N 3 ) 
1 

IP' 



2n + )a 



N 3 )bc--^a(b 2 + c 2 



+{q - 1 - AK + 4£^ 
_1(7V 1+ Ar 2 
b" = (q-2)b' 

+(q-l-AK+ (E+ + V3E_)(E+ + V3E_ + 2) - N 2 (N 3 + Ni) -11+ - V3U-)b 
~(Ni + N 2 + N 3 )ca - j^b{c 2 + a 2 ), 
c" = (q-2)c' 

+(q-l-4K+ (E+ - v/3E_)(E+ - V3E- + 2) - AT 3 (iVi + N 2 ) -11+ + V3U-)c 
-L(N 1 + N 2 + N 3 )ab-j^c{a 2 + b 2 ) . 



(2.26) 



(2.27) 



(2.28) 



The basic variables are therefore [H, E + , £_, Nx, N 2 , N 3 , a, b, c], and they are determined uniquely and 
completely by Eqs. (E~H|)-(EIIHJ), (|2~2^|) - (|2~2^|) . along with the constraint equation (|?TT2"|) . 
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III. BIANCHI I CO-EVOLUTION OF YANG-MILLS FIELDS AND PERFECT FLUIDS 



We compute the solutions of the dynamical equations (|2.14|l - ll2.19ll and l|2.26[l - l|2.28[l numerically and illustrate the 
typical evolutionary behaviours. In this section, we consider for simplicity only the Bianchi type I spacetime (so 
if = 0) with black-body radiation (7 = 4/3) or dust (7 = 1) and a YM field. We set initial value of H as H in [ = 10 -10 
(our unit), which is not a restriction because our time unit depends on <?ym which determines the physical identity of 
the YM field and can be scaled. 

First, we display the behavior of the three leading quantities, S 2 , Oym, and il m . Figure shows S 2 — > 0, i.e., 




FIG. 1: The typical behaviors of E 2 , Hym, and f2 m with increasing time. This figure shows how E 2 — > 0, and the spacetime 
approaches the isotropic Friedmann model. It is remarkable that JIym does not damp. 



spacetime approaches to the fiat isotropic FRW universe. We also find that SIym does not damp, but approaches a 
constant. It is interesting to compare this behavior with that which obtains close to isotropy in universes containing 
radiation and a magnetic field. There we found that £ and f2 ma g both fell slowly as (lnt)^ 1 as t increased, with 
asymptotic approach to an attractor where Slmag/^rad = AS, where A ~ O(l) is a calculable constant depending on 
the orientation of the field and the number of effectively massless spin states in the radiation background. 

Second, we can display the evolution of the components of the YM field. Figure [3 shows that a, 6, and c follow a 




FIG 2: Typical behaviors of the Yang-Mills field. 

(a) a, b, and c damp whilst undergoing chaotic oscillations of increasing frequency. 

(b) The averaged damping rate of a, b, and c has the approximate form e~ T . 



sequence of damped chaotic oscillations, and decay as e _T . Figure|21also shows that the frequencies of the oscillations 
of a, 6, and c are growing, although their amplitudes are decaying. This is further supported by Figure |3 which 
shows that the mean amplitudes of a', 6', and c' remain nearly constant. 
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FIG. 3: Typical behaviors of the time derivative of the principal components of the Yang-Mills field. The mean amplitudes of 
this shows amplitudes of a' , b' , and c' remain nearly constant. 



Third, we show the behavior of II-|-, the Hubble-normalized anisotropic pressures. Figure^shows that the Tl± oscil- 

0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 

012345678 
x 

FIG. 4: The typical behaviors of the anisotropic pressure components II±. This shows that II ± oscillate in a complicated way 
with nearly constant mean amplitudes. 

late in a complicated fashion with approximately constant mean amplitudes. This is consistent with the observation 
that Oym cx fi m but it is surprising that the pressure anisotropy remains constant despite the fall in the shear. These 
chaotic oscillations are a familiar feature of the YM evolution which resembles the chaotic Hamiltonian dynamics 
of a point bouncing inside a potential whose four steep walls in the x — y plane are formed by the branches of the 
rectangular hyperbolae x 2 y 2 — constant [TcHll). 

In the next section we try to understand these features emerging from the numerical studies of Bianchi I spacetime 
with the YM field and radiation or dust. 




A. The approach to isotropy 

In this section we consider Bianchi I spacetime again, i.e., N a — and K = 0. We assume that the FRW limit 
holds, that is £ — > 0, so the dynamics can be close to isotropy. This results in 

9 = Y M + ^(37-2)O ra , (3.1) 

S^YM + f^m = 1 ■ (3-2) 
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These conditions imply 



q= l+^(3 7 -4)fi m . (3.3) 



Motivated by the results of the numerical analysis, we also assume that £Iym ~ > const, and f2 m — > const. These 
assumptions reduce <|2.25[1 to the condition that 

n' m = {2 g -(3 7 -2)}Q m = 0, (3.4) 

and if Vt m ^ 0, we have 

g=i(37-2). (3.5) 
Hence the evolution of is given by if' = — (1 + q)H = — § 7 if. Solving this equation, we have 

H cx exp ^— ~jl 7 ^j ■ (3-6) 

This means that to leading order the YM field does not determine the expansion rate in FRW limit. 
From eqs. I|3.3[l and (|3.5|l . we get 

3 7 - 4 = (3 7 - 4)O m . (3.7) 

Therefore, in the radiation case ( 7 = 4/3), f2 m can converge to any value. Otherwise ( 7 ^ 4/3), f2 m becomes dominant, 
and we must have Q m — ► 1. In the radiation case ( 7 = 4/3), the results become 



q = I, 0,ym + = Ij H oc e 



-2t 



Hereafter, we will consider only the radiation case. 

From our numerical results, we know that the YM field oscillates with growing frequency while its amplitude is 
damped. Therefore, we assume that a — > 0, a' -f* 0. Because the U± are both finite 3 , eq. (|2.26ll becomes 

a" = -a' - -^a(b 2 + c 2 ) = -a 1 - e 4T a{b 2 + c 2 ) . (3.8) 

If we divide a into an amplitude and a phase, with a = e~ pT e 1,ei , where p is the damping rate and 9\{t) is the 
time-dependent oscillation frequency, then 

a' = -pa + i0[a = {-p + i0[)a (3.9) 
In order to keep a' 0, 9[ has to grow like 

8[ ~ 0{e pT ) . 

Differentiating eq. I|3.9|l . this gives 

a" = -pa' + i9'{a + i9' 1 a' 

= {p 2 -0' 2 +i{0'l-2pe l 1 )}a, (3.10) 

where 9[ 2 grows 0(e 2pT ) and 9'{ depends on the phase of 9[. Substituting in eq. 1)3. 8|) yields 

p 2 - 9' 2 + i{9'l - P 8' l ) = -{-p + i9[) - e 4r (b 2 + c 2 ) 

= p - i9[ - e iT e- 2pT {e 2t92 + e 2 ^ 3 ) , (3.11) 

where, for the simplicity, we assume b and c damp at the same rate as a, so b = e~ pT e t92 and c = e~ pT e l0:i . Picking 
out the dominant terms in eq. I|3.11[) . we have 

-9' 2 +l9'l = _ e (4-2p)r (e 2tf 2 + e 2ie 3) 

= -e (4 - 2p ^ T 2cos(8 2 ~8 3 )e^ e2+e;l) . (3.12) 



In Bianchi I, II, and VIo spacetimes, —1 < il+ < 2 and — %/3 < II— < v'S (see llll t. 
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It is generic to assume that there is no relation between 62 and 6*3, so 62 63. As a result, the characteristic frequencies 
of the above equation will be neither purely real nor purely imaginary. Therefore it is likely that 0(6[ 2 ) ~ 0(6"), 
which means that the leading-order solution of the above equation grows as e 2pT , i.e., p=l. 



From this, we can estimate H± as follows, 



n. 



1 



= -! ( - 2 a' 2 + b' 2 + c! 2 - 2e 2l ^ +e ^ T + e 2i ^ +e ^ T + e 
6 



2i(e 1 +e 2 )T 



O (const.) 



IL 



O (const.) 



(3.13) 



V3 / 
6 V 

It follows from our discussion that II± y4 despite the fact that E — > . This result is quite unexpected. 

It is worth checking explicitly that these non-zero anisotropic pressures, H±, do not break isotropy. That is, we 
need to confirm whether these non-zero II ± make S 7^ or not. Using eqs. (|2.23|l and l|2.24|l . and the property that 
E± < n±, and assuming A 2 + A 2 ~ B 2 + B 2 ~ C 2 + C 2 , where A = a' + (-2E+ + l)a, B = b' + (E+ + V3E_ + 1)6, 



C 



(Eh 



V3E_ 
E' 



- l)c, A 
(E 2 )' 



iV lfl + f , B = iV 2 6 + 



1 



2E 



(E + E^ + E_E'_) 



1 

E 



(A 2 + A 2 ) cos * + (B 2 + B 2 ) cos ( * 
(.4 2 +yl 2 ) I cos* f cos j »!' 



and C = iV 3 c + 

(E+n+ + E_n_) 

(C 2 



w, we have 



2 

— 7T 

3 



T 2 

* 7T 

3 



0. 



(3.14) 



where \J> = tan 1 (E_/E + ) . This shows that the non-zero IIh- modes do not break isotropy. The assumption, _4 2 +.4 2 
B 2 + B 2 ~ C 2 +C 2 , is justified by numerical solutions. Figure [S] shows the time averages of A 2 +A 2 
over a period of O.lr, which confirms our assumption. 



B 2 +B 2 ~C 2 +C 2 



0.5 




FIG. 5: The typical behaviors of the (A 2 + A 2 ), (B 2 + B 2 ), and (C 2 + C 2 ), where (• ■ ■) means time average per O.lr. This 
shows that A 2 + A 2 , B 2 + B 2 , and C 2 + C 2 are nearly equal after the shear is significantly reduced, i.e. r > 2, although 
A 2 , B 2 , C 2 , A 2 , B 2 , and C 2 oscillate strongly. 



In summary: so far our analyses have shown that if we assume that the radiation- YM universe approaches FRW, 
in the sense that E — > 0, and also that Hym ~~ ¥ const, and Q m — > const., then the YM field components satisfy a — > 0, 
and a' -> (9(const.), with A 2 + A 2 ~ B 2 + B 2 ~ C 2 + C 2 , and it is found that H cx e _2T , a cx 0(e~ T ), E' -> 0, and 
the normalised pressures II± oscillate with nearly constant amplitudes. In the next section we provide some further 
details of this evolution which can be obtained from the numerical studies. 



IV. A SURVEY OF TYPICAL EVOLUTIONARY BEHAVIORS 



Our numerical studies of the Bianchi I expansion dynamics were performed for various combinations of radiation 
or dust and YM field. We give here a summary of the principal conclusions from these investigations. These results 
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are stable against small changes in the initial data (-ffini and Ej n i) and appear to be robust. This robustness is a 
consequence of the generic nature of the chaotic oscillations of the YM field and the fact that these chaotic oscillations 
have a small effect on the expansion dynamics of the universe. 



A. Radiation and YM field 



1. YM field dominates: Uym 3> f2 ra d 

In the case where the YM field density dominates the density of the isotropic black-body radiation field, we find 
that Oym ~ const, and O rac j ~ const, during the evolution. If we define an averaged expansion scale factor for the 
expanding universe, £{t), by 

TT _ dijdt 

then the normalised shear is found to fall as E oc £~ 1 - 01 in the numerical integrations, while the Hubble expansion 
rate falls at the rate expected in a FRW model, with H oc £~ 2 - 00 , so the mean scale factor evolves as I oc t°- 500 ~ i 1 / 2 j 
as in a FRW universe. 



2. Radiation dominates: Qym <C firad 



If the black-body radiation density dominates the YM field then we still find Uym ~ const, and f2 ra d ~ const, and 
the normalised shear falls almost linearly with the scale factor, as E oc ^ _107 ,with H oc £~ 2 - 00 and £ oc £ - 500 ~ i 1 / 2 as 
in the FRW model. Thus we see that the evolution in both of the radiation plus YM field situations is similar, with 

E oc I' 1 and t oc i 1/2 

holding to an excellent approximation. Although the anisotropic pressure can be significant, it is not driving significant 
shear expansion anisotropy. 



B. Dust and YM field 



1. YM field dominates: Qym S> ^Wst 

When the YM field dominates over the dust we are in the very early stages of the overall evolution because on 
average the dust density redshifts away more slowly than the YM field. When the YM field still dominates, so 
Oym ~ 1, the numerical evolution gives 

ttdust OC ^°- 983 

e oc r 1 - 01 

while H oc £~ 2m so £ oc t ' 500 ~ t 1 ^ 2 . This reflects the assumption of the YM field domination and is not inconsistent 
with eq. 1)3.6(1 because the assumption used there, IIym and flm — > const, is violated. We see that the damping 
rate of shear and the expansion law are the same as that in the radiation-dominated case. It is understandable that 
^dust oc ^°- 983 , because /idust oc £~~ 3 and H oc t^ 1 oc £~ 2 . This results in the dependence rWst = Mdust/3-ff 2 oc £, 
as expected. Thus, f^dust grows and will eventually become dominant and the assumption that J]ym 2> ^dust will 
eventually fail. 



2. Dust dominates: Qym <K ^Wst 



This is the natural situation for the universe to evolve into after the radiation-dominated era. Our numerical studies 
find that in the dust-dominated phase the YM density falls as 

o ~- H-0.954 
SiYM OC I 
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while the normalised shear falls as 

e oc r 1 - 48 ~ r 3 ' 2 

and the Hubble rate falls as if oc £~ 1 - 52 . Hence, we obtain a close approximation to the evolution expected in a 
dust-dominated FRW universe, with I cx £ - 658 ~ £ 2 / 3 . As expected, the shear falls off more rapidly than in the 
radiation-dominated case because of the growing influence of the isotropising dust density. It is also interesting that 
the shear falls off more rapidly than in the simpler dust plus magnetic universes Q . We note that the fall of the shear 
in the YM case is close to the E oc £~ 3 / 2 fall-off that would occur if the YM field were absent in a dust-dominated 
Bianchi I universe. We see how this arises by noting that 

EV = -^E + +H + ~-^E + 

because |n±| < Qym, and therefore, 




V. COMPARISON OF YANG-MILLS AND MAGNETIC FIELDS 

One of the original motivations for our study of the evolution of YM fields in the presence of perfect fluids was 
the unusual behavior found in the case of a pure magnetic field and a perfect fluid, notably in the situation where 
the perfect fluid is black-body radiation. Since the YM field is a generalisation of a magnetic field it is instructive to 
compare and contrast the results for these two cases. 

In the case of an almost isotropic universe (I oc t 1 ^ 2 ) containing a pure magnetic field (or other anisotropic stresses 
with pressure anisotropy of the form l|l.l|0 and black-body radiation with il mag <C fi ra d, the expansion- normalised 
shear falls logarithmically in time, and 

E OC /imag/Mrad OC 1/ In t OC 1/ 111 £ 

This unusual 'critical' evolution arises because there is a zero eigenvalue when we linearise the Bianchi type I magnetic 
radiation universe around the isotropic Friedmann radiation universe 0, Q . Note also that the ratio of the magnetic 
to the black-body density is not constant as is often assumed, but falls slowly due to the coupling to the shear 
There is a late-time attractor with /i mag //irad ~ 5] and 

^mag = Mmag/3-ff" OC /^mag/Mrad OC E. 

It is interesting that f2 mag oc E. This means t licit ^^niag 

can be constrained directly by the CMB temperature 
anisotropy, AT/T, since the presently observed AT/T ~ E rcc where z ICC ~ 1100 is the recombination redshift. After 
accounting for the short period of dust dominated evolution from the equal-density redshift, z cq ~ 10 4 , to z roc this 
leads to strong limits on any spatially homogeneous cosmological magnetic field today of 3.4 x 10 _9 (rio^5o) 1 ^ 2 G, 0, 
or on any anisotropic stress with pressures of the form (|l.lf) . see ref. [4J for details and examples. 

In a dust-dominated era (f2 mag <C ^dust) the magnetic stresses still slow the fall off of the shear to E oc l~ x oc i~ 2 / 3 , 
whereas we would have had E oc £~ 3 / 2 if the magnetic field was absent. The magnetic density evolves as 

^mag OC /^mag//^dust OC £ 

The results for the different magnetic and YM evolutions are summarised in Table [I] 



Material content 


Shear Evolution, E 


Anisotropic stress density 


Dust + magnetic field 


E oc r 1 oc t-' 2/s 


Sitnag S CX t ^ 


Dust + YM field 


E oc £~ 3/2 oc t- 1 


n Y yi oc t~ 2/3 


Radiation + magnetic field 


E oc (ln^)- 1 cc (hit)- 1 


fi„ H g oc E oc (hit) -1 


Radiation + YM field 


Eoc^" 1 oct" 1/2 


Aym = const. 



TABLE I: Evolution of the normalised shear, E = a/H. 
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VI. OBSERVATIONAL BOUNDS 
A. The microwave background 

The key difference that these results display between the magnetic and YM fields is the rapid fall-off in the shear 
that accompanies the YM evolution. As a result the YM field density does not determine a shear attractor at the 
end of the radiation era and does not have a significant effect on the CMB anisotropy in the way that the magnetic 
field does. As a corollary, observations of the CMB temperature anisotropy do not provide a direct and powerful 
upper bound on the present YM field density in the way that they do for a magnetic field density. A scenario with 
which /iYM ~ A'rad is not constrained by shear anisotropy of CMB, although it may be constrained by Big Bang 
nucleosynthesis (BBN), as we discuss below. 

The slow decay of the shear in the magnetic universe means that the CMB anisotropy places a stronger limit on 
the magnetic field density than can be obtained by considering its effects on the expansion rate of the universe at the 
time of BBN, zbbn ~ 10 9 — 10 10 . However, in the YM case the rapid fall in the shear as we go forward in time means 
that a negligible shear at the the time of recombination will be much larger at earlier epochs, with an enhancement 
factor of 

Ebbn ^ 1 + zbbn ^ ^5 

^rcc 1 ~\~ Zeq 

Since E BB n < 0.35 from the helium-4 abundance pi IM I2I I23] . we see that the BBN anisotropy bound is 

S rcc < 3.5 x 10~ 6 . (6.1) 

This is about an order of magnitude stronger than the limit imposed by the CMB since Q 

AT/T = £ rec x f(9,4>,Q ), (6.2) 

where / ~ 0(1) is a dimcnsionless pattern orientation factor, which is bounded as 0.6 < / < 2.2 in flat or open 
universe by the COBE results over > 10° angular scales and so the Bianchi I CMB constraint on E rcc is only 

E rec < 10™ 5 . (6.3) 

As discussed in Barrow |24| . this simple anisotropy evolution leads to an unphysical super-Planck anisotropy energy 
density in the shear modes at very early times if extrapolated backwards to z 3> z cq from modest values of shear at 
z oq . A physically reasonable requirement (the 'Planck Equipartition Principle (PEP)' would be that at t p i ~ G 1 / 2 ~ 
10 _43 s all energy densities contributing to the cosmological dynamics should be bounded above by the Planck density 
fipl ~ t~j ~ 10 94 <?m.cm~ 3 . This would require E p / < 1 at z p i ~ 10 32 and hence 



1 + z r 



-'pi 1 + z pi 
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and the residual anisotropy in the CMB on large angular scales in the YM Bianchi I universe would be completely 
negligible. By contrast, in the magnetic universe the slow logarithmic fall in £ allows an interesting anisotropy level 
Srec ~ O(10 ) to persist at recombination even if E p ; < 1. The predicted value of AT/T ~ E roc depends upon 
the number of relativistic spin states contributing to the equilibrium radiation sea and can also differ if other forms 
of anisotropy are assumed 24] . In general, we would expect that YM fields in the more complex Bianchi type VII 
universes would display slow logarithmic decay of their shear regardless of the presence or absence of the YM fields. 
This is due to the anisotropic 3-curvature which mimics the presence of an anisotropic 'fluid' of l ong - wavelength 
gravitational waves satisfying (|l.l|l to a good approximation and leads to E oc (lnt) -1 , see refs. [25l l26l l27t p g l [29 | . 
Note however, that these modes are excluded in anisotropic open universes if their spatial topology is compact |3(ll3l|: 
hyperbolic Bianchi spaces with compact topology must be isotropic. 

B. Constraints from BBN 

In the scenario containing YM fields studied above, the YM density Oym does not evolve in proportional to the 
shear. The YM pressure anisotropy does not dominate the simple 'adiabatic' decay of the shear (a cx £~ 3 ) that occurs 
in the absence of anisotropic sources or anisotropic 3-curvature. Therefore it is impossible to constrain Oym by the 
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shear anisotropy alone, in the way that f2 mag was constrained by shear. Primordial nucleosynthesis gives the strongest 
constraint on Qym- I n our YM scenario, the expansion rate H is larger than in the isotropic FRW case, because 

H 2 = (8vrG/3)(AWi + /iYM) 

where /i ra( j oc T 4 is fixed by observation. The increase in the expansion rate raises the temperature at which neutron- 
proton abundance ratio freezes out of equilibrium, and this leads to an enhancement in the final helium-4 abundance. 
The YM field evolves on average like a radiation field and we can use constraints on the density of dark radiation 
(DR) 32] to limit its allowed effect in this process. This gave a bound of S1 DRbbn < 0.105fi ra d BBN at 2a confidence 
level. Therefore, we expect the YM field density to be similarly constrained with the bound 

Hym < 0.105f2 rac i. 

A similar bound would be obtained from BBN considerations for the magnetic energy density. However, the bound 
previously obtained from the CMB anisotropy is far stronger in this case: f2 mag < 10 -5 r2 ra d- 

VII. CONCLUSIONS 

We have analysed the evolution of YM fields in anisotropic radiation and dust cosmologies which evolve close to 
isotropy. The YM field undergoes chaotic oscillations which produce small chaotic vibrations about the Friedmann 
expansion when the anisotropy level is small. We investigated the evolution of the shear anisotropy and the YM 
density during the radiation and dust eras. This revealed significant differences to the unusual situation that is known 
to exist in magnetic cosmologies containing perfect fluids. In particular, unlike in magnetic universes, there is no 
attractor in the radiation era which couple the YM energy density to the shear anisotropy. Consequently, there is 
no direct bound on the YM field density from the CMB temperature anisotropy. We have carried out a comparative 
analysis of the magnetic and YM universes which shows how magnetic universes of Bianchi type I are principally 
constrained by the CMB anisotropy whilst the YM universes of this type are constrained by the effects of the YM 
energy density on the synthesis of helium-4. The YM evolution reveals unusual features. Despite the fall in shear to 
Hubble expansion ratio, the Hubble-normalised pressure anisotropies induced by the YM field do not decay during 
the radiation era, but oscillate chaotically with constant amplitudes. 
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